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Abstract 

The derivative expansion of the one-loop effective Lagrangian in QED4 is con- 
sidered. The first term in such an expansion is the famous Schwinger result for 
a constant electromagnetic field. In this paper we give an explicit expression 
for the next term containing two derivatives of the field strength F^,^. The re- 
sults are presented for both fermion and scalar electrodynamics. Some possible 
applications of an inhomogeneous external field are pointed out. 
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1 Introduction 



The effective Lagrangian for a constant electromagnetic field was calculated by 
Schwinger more than forty years ago []T|. The Schwinger Lagrangian has the virtue 
of being nonpertubative, taking into account the effects of the background field 
to all orders. However, its validity is restricted to constant fields and, since that 
time, much efforts have been devoted to finding a systematic way to generalize 
the result for strong slowly varying fields P, ^ or for the fields localized in a 
finite region of space 0]. Such a generalization would be of great interest from 
the physical point of view, since it can be applied to studying instabilities (such 
as spontaneous symmetry breakdown) in quantum field theory in external field 
or it can be used for the explanation of the narrow e~e^ peaks observed 
in heavy-ion scattering experiments ||^, where strong varying electromagnetic 
fields are indispensable ingredients. Also, nonperturbative information on the 
fermion determinant, closely connected with the one-loop effective action, is very 
important because fermion determinants produce an effective measure for gauge 
fields when the fermions are integrated out p]. 

In the present paper, we calculate the low-energy effective action for fermions and 
bosons in QED4 as a derivative expansion around nonvanishing field strength 

-^e// = ^q{Fhv) + d\Fal3d.yF„si^l"^^"^ {F^u) + • • • , (1) 

where £fc's are some local functions of the field-strength tensor F^^. In the 
derivative expansion (|^) the first term Cq is the famous Schwinger result and we 
compute the next term Ci. 

Recently, an analogous calculation has been performed for QED in 2 + 1 dimen- 
sions . As for 3 + 1 dimensions, we mention the paper [Q] where the heat-kernel 
method was adopted to obtain the derivative expansion (p. However, it is impos- 
sible to compare our results with those in because some important functions 
in the cited paper were not given explicitly. We comment more on comparison 
with results of Refs. [0, ^ at the end of the paper. 
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2 Derivative Expansion of the One— Loop Effec- 
tive Action in QED4 



The one-loop effective action in QED reduces to computing the fermion deter- 
minant 

ziy (1) (A) = In Det(zP - m) = i In Det (^-Pj + ^(^i.uF'"' + 

= i 1 d^xixltr In (^-Vl + |cT^.i^^" + \x). (2) 

Here V = 'y'^V^ and = + ieA^ is the covariant derivative, a^^, = i/2[y^, ■ji,] 
and tr refers to spinor indices of Dirac matrices 7^; \x) are eigenstates for some 
self-conjugated operators x^ of coordinates. Throughout the paper we use the 
Minkowski metric rjn^, = (1, —1, —1, —1). 

With use of the identity lnA = — /q°° exp{—iTA)dT/T for introducing the proper- 
time coordinate, Eq.(|^) can be represented through the diagonal matrix elements 
of the operator U{t) = exp(— zrif): 



00 



lyW(A) ='-J d^x J ^e-''^'^tr{x\exp{-tTH)\x), (3) 


where the second order differential operator H is given by 

H = U,U'' + ^a,,F>^''{x), U^ = -iV,. (4) 
The evolution operator U (r) satisfies the Schrodinger equation 

= HUir) (5) 

together with the initial condition t/(0) = 1. In Euclidean space this operator is 
closely connected with the heat kernel operator. Because the matrix elements of 
U (r) cannot be calculated for arbitrary external fields, several approximate regu- 
lar schemes have been developed: the technique of asymptotic expansion at small 
values of the proper time (r — 0) [9-13], covariant perturbative theory [|I^ |15 



and, at last, pseudodifferential operator technique [16-19]. The asymptotic ex- 
pansion at r — i> is equivalent to the asymptotic expansion of the effective action 
in inverse powers of mass parameter m^. The coefficients of this expansion are 
the space-time integrals of local invariants, -E„(x, x) {n = 0, 1,2, . . .), of grow- 
ing powers in field strenth F^,^ and its derivatives, called DeWitt-Seeley-Gilkey 
coefficients. The local (Schwinger-DeWitt) expansion is, therefore, an asymp- 
totic approximation valid for weak and slowly varying background fields. The 
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covariant perturbation theory of Barvinsky and Vilkovisky deals with the case of 
weak and rapidly varying fields. It corresponds to summing up all terms in the 
Schwinger-DeWitt series with a given power of field strength and any number 
of its derivatives. The covariant perturbation theory like the small proper-time 
expansion, fails in the case of large fields. An approximation scheme complimen- 
tary to the above mentioned ones would correspond to strong and slowly varying 
fields (derivative expansion). The lowest-order approximation of such a scheme is 
the case of covariantly constant background fields. In principle, the most general 
approach of pseudodifferential operators permits one to develop such a derivative 



expansion (see, for example, and |@]), but in the present paper we shall follow 



another approach to derive the expansion (|T]). 

The matrix elements {x\exp{—iTH)\x) entering the right-hand side of Eq.(|^) 
may be treated as the matrix elements of the evolution operator of a spinning 
particle, r and H being the proper time and Hamiltonian of the particle. The 
corresponding canonical momenta are P^'s which obey the commutation relations 
[x^,P'^] = i5^^ and (x|P^|?/) = —id^6{x — y). Following the standard method ^1 



the transition amplitude {z\U{T)\y) between points x(0) = y and x{t) = z can be 
written in terms of a path integral over c- number (x^(t)) and Grassmann {'ip^{t)) 
coordinates: 

tr{z\U{T)\y) = N-^ f V[x{t),^{t)]exp U j dt[Lbos{x{t)) + Lf,r{Ht),x{t))] 





(6) 

where is a normalization factor, and 

r / \ 1 dx doc A / \ dx / 

Lf,r{i^, x) = \'^u'^- ien^^'Fuxix). (8) 

The integration in Eq.(^ goes over trajectories a;^(t) and ^/'^(t) parametrized by 
t G [0, r] and obeying the boudary conditions 

xiO) = y, xiT) = z, ^(0) = -^(r). (9) 

The Lagrangian L^osix) + Lfer{'ip,x) may be treated as Lagrangian of the rel- 
ativistic spinning particle in external electromagnetic field The 
first-quantized representation for the effective action, given by Eqs.(^, (|)-(||), 

^For a more sophisticated representation in terms of path-integral with reparametrization- 
invariant Lagrangian, see ||2j, We note that from the point of view of further using the first 
quantized path integral in applications it would be quite desirable to have a representation for 
the evolution operator in terms of boson operators only. Since 7-matrices enter the operator 
(0) through cr^y, satisfying 5*0(3, 1) Lie algebra, such a representation should exist. 
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can be used as a starting point for developing approximate schemes for computing 
the expansion in the case of arbitrary external fields. The path integral (^) is the 
particle theory analogue of the Polyakov path integral in string theory, and mo- 
tivated by this fact, Strassler |^ has developed the first-quantized perturbation 
theory for calculating the one-loop effective action. Moreover, such an approach 
to ordinary field theory, based on path integrals defined on one-dimensional world 
lines, can be extended to the evaluation of Feynman diagrams for Green's func- 
tions at arbitrary loop order ||25| and to the computation of form-factors p6 |. 
In the last case it was shown how the higher derivative terms appearing in the 
effective action can be effectively summed up and lead to the results obtained in 



the heat kernel approach of Barvinsky and Vilkovisky [14 



In order to generate the expansion (|1|) from the path integral it is convenient 
to choose a special gauge condition for the potential A^{x), namely the Fock- 
Schwinger gauge 



(10) 



which leads to the series 
1 



AJx) 



{x' - y^)FUy) + ^{x^ - y^){x'' - ynd.Fxu{y) 



E 

n=0 



(x^ - y^){x''' - y"^) . . . {x"" - y"") 
n\{n + 2) 



du,du^...d^^Fxu{y). (11) 



Carring out the change of the variable x{t) for x'(t) = x{t) —y in the path integral 
(I) (henceforth we omit the prime) and substituting (|lT]) into (|^), we obtain 



tr{z\U{T)\y) = N~ 



X— + L,(x) 



exp 



D[x{t),ij{t)]exp 
dip 



, / 1 dxy dx^ 
dt\ — 



4 dt dt 



I 

^ 



x'-FUy) 



dt ( 2^-^ - ^erVFux{y) + L^{x, ^) 



(12) 



The new boundary conditions for x{t) are x(0) = and x(r) = z — y. In Eq.(^) 
F^i, does not depend on x{t). The expressions for the interacting terms, L2{x) 
and L^^x^ip), containing derivatives of F^i, with respect to coordinates take the 
form: 



L2{x) 



dt 



°° eF dr^o 
ti n\{n + 2) 
e dx''' 
3 "''''11 



x"' (t)... a;^"+i (t) 



x^x" + ^F, 



dx^ 



u\,aK ' 



dt 



■x^x^x'^ + ... 



(13) 
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n=l 



le 



Here we use a conventional notation for derivatives: F: 



du^Fx^{x). 



Introducing c-number and Grassmann external sources, we rewrite 



i j dtlL2 



1 6 



i6ri{t) 



tr{z\U{T)\y) = exp 
where the generating functional is defined by 



1 6 
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(14) 



as 



(15) 



ZAv,^]{z;y) = / D[x{t),tlj{t)]exp 







1 (^iJy jy (^tXy 

2 dt dt 



- ex^Fxuiy) 



dx" 

X — + 2ri,x' 
dt 



exp 



^ ^ 



(16) 



Since the expression in the exponent of the right-hand side of (|T6[) is quadratic in 
the variables x and if), the calculation of the generating functional Zr{ri,S,]{z]y) 
can be reduced to the calculation of the determinants of the following one- 
dimensional differential operators: 



^ 2 e^A^-^^ 



and O2 = ^VfMi^^ - '^ieF^u, 



(17) 



defined on the interval [0, r] with the periodic and antiperiodic boudary condi- 
tions, respectively. Taking into account the initial condition 



(18) 



which is equivalent to the operator equality U{0) = 1, and introducing two inde- 
pendent invariants constructed from the field tensor F^j^: 



and 



with 



-^F^'-'F^, and g = ^e^-^^^F^^F, 



(19) 



(20) 
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one can represent the result of the integration in (|16|) at coincident arguments 

z = y = X a.s 

ZtYi,^[x]x) = —-———-{eTK^){eTK_)coth.{eTK^)cot{eTK_) 



■ exp(^^,fM-l5r[e]), (21) 
where the quadratic in external sourses forms, S^f^ and S'Jf^, are defined by 

r r 

Scriv] = JdhJdt2Vu{ti)D''\h,t2)vx{t2), (22) 


r r 

S!r[^] = JdhJdhUti)S''\h,h)Ut2), (23) 



The Green functions entering the last expression are defined by the formulae: 
+ coth(e/jr) (1 + exp[2e/,(ti - ts)]) 

- sinh-i(e/,r)(exp[e/,(r-2t2)]+exp[e/,(2ti-r)])j, (24) 
S^'ih, h) = ^ ^0)1 - t2) - tanh(e/,r)) exp[2e/,(ti - ts)]. (25) 



Here we use the matrices A'^^^ which were originaly introduced in [^]. The explicit 
expressions for them are 

- ^^J^^ , (26) 



where 



The main property of matrices (p6D that was used in (pi]) and (]25|) is 



F''^A^i)x^^ = = fiA'^i)^. (27) 

Other useful properties of these matrices that will be used below are: 

E 4) = <(0 = 1' 4)^^A0-) = S,jA%y (28) 

i=l 



As can be easily seen, in the case of vanishing field the propag ators D''^{ti,t2) 



and S{ti,t2) are reduced to the free Green functions used in p6 



Substitution of ( PT] ) and (P^, (^) into (|T3p leads to the expression for tr{x\U\x). 
The last step is to expand the exponent in terms of derivatives and to keep only 
two-derivative terms. Therefore, one obtains: 



tr{x\U{T)\x) 



, T 1 

1+^1 f/t[V2(t) + W2(t)] - 2 / / ^^l^i2[Vi(ti)Vi(t2) 
' 



+Wi(ti)Wi(t2)] - // c/tirft2Vi(ti)Wi(t2))^.[r7,e](^,a:) 



where 



Vi(t) 
V2(t) 
Wi(t) 

W2(t) 







— \ lim 

-(^Fyx,i,K lim 



(29) 



,7=0,^=0 



- F ^ 



(30) 



Substituting the generating function ( pT]) which depends on the Green functions 



(|J) and (p5D , one can rewrite (^) in the form: 

i{eTK^){eTK^) cot(erfC_) coth(eri^'4 



tr(x|t/(r)|x) 



1 g ^Fy\,^K 



T 

X y dt(^D''\t, t)D''^{t, t) + I)'^^(t, t)D^''(t, t) + I)'''^(t, t)D^^(t, t) + 45'^^(t, t) 


T T 



-25"\2, 1)S'"'{2, 1)) + 65"''(2, 2) (l)''^!, 1)D^^(1, 2) + Z)'^^(l, 1)D^^(1, 2) 
+I)'^^(1, 2)D^^{1, 1)) + I)'^^!, 1)Z)"''(2, 2)D^^(1, 2) + 2Z)'^\l, 1) (l)"^(2, 2) 
X D^''(l, 2) + Z)'"'^(2, 1)L)'^^(2, 2)) + D'^'^il, 1)D''p{2, 2)D^''{1, 2) + 2I)'^'^(1, 2) 
X (l)"''(2, 2)D^'^(1, 1) + D'"'{2, 1)D^^{1, 2)) + D'^^'il, 2)Z)"\2, 1)D^''(1, 2) 
+D''P{1, 2)Z)"'^(2, 1)D^"(1, 2) + D'^'il, 2) (d^'^(1, 1)D'^^(2, 2) + D^^{1, 2) 
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Here the dotted functions are defined by 



(31) 



D^%1,2) =^ (32) 

D^^(l,2) _^Z}M^(t^,t2), (33) 

Z)^^(t,t) \ir^J-D^^(t,,t). (34) 

The straightforward, though tedious, computation gives the result 
tr{x\U{T)\x) = - , I ^ {eTK_){eTK+) cot(erir_) coth(erir+) 



^ 3,1 ^ 

-C^2^''(/,, fu fk) Ko^(o4') + ^S4)4^) + 2^(i)4)4'); 



(fc) 



(35) 



where the coefficients C^^ {a, (where X, Y E {V,W}) are given in Ap- 
pendix. 



In principle, the expression (|35| ) is our final result. It remains only to substitute 
the last expression into (|^) and to renormalize field and charge according to the 
rule 



{^'"'f = YTc^'' = v/T+C^^M, (36) 



oo 



C = ^/^exp(-zmV). (37) 







Note that the renormalization used merely means performing some subtraction 
(precisely the same as in the origional Schwinger paper |^) from the term contain- 
ing no derivatives of field strength and changing all bare quantities for physical 
ones. 
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As was mentioned at the begining of the paper, the derivative expansion in QED4 
was also studied in 0]. The result of that paper, however, was presented in an 
explicit form only for the special case of the electromagnetic field: 

g = 0, F'^'ix) = <l>(x)F^^ (38) 

where $(x) is a scalar function and F^'^ is a constant matrix. Under the particular 
choice of the field configuration given by (^) the derivative part of the diagonal 
matrix elements (^) reads 

tr{x\U{T)\x)aer = "Ji (3^^^ " ^coY' " Acu^Y^ + 3coY + cu'), (39) 

where u = erV2JF, Y = coth(ci;). We note that this result was obtained in by 
means of symbolic calculations on a computer using the REDUCE package. 

Now turning to the calculation of the derivative expansion for the scalar electro- 
dynamics, one does not need to repeat all the calculations similar to those made 
above. In order to see this, we recall that the effective one-loop Lagrangian in 
this case reads 



oo 



42(x) = ~^ I -{x\U,Ur)\x)e-^"^'\ (40) 







The evolution connected with the transition amplitude, {z\Ubos{T)\y), is described 
now by the Hamiltonian (compare with (^ and (^)): 

Hbos = {^^^ - eA^{x)){td^ - eA'^(x)). (41) 

Thus, omitting all terms originating from the fermion part in (^, i.e. putting 
L3 = in ([T^), (p!5| ) and S^l'^ = in (|21|), we come to the following expression 

{x\Ubos{r)\x} ^ 



IGvr^r^ sm(eTK-) sinh(eri^'4 



X 



X 



'^^rf. f, f,']( A''^ A'^'^ Af^'^ -\- A'^^^ A'^f A^'^ -\- 9 A'^^ A'^P A^^'^ 

+2cr{f,, fi, fk) (^S4)45 + ^^uAiAk)) + c^'^if^^ f^)^u)4)^ik) 



(42) 
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The coefficients used here are the same as those that were used in (|35|) . In the 
special case of the electromagnetic field of the form (|38|) we obtain 



/ / M V 1 (9.$(9'^$ /3(cucothcj- 1) , \ uo 

{x\Uhos[r)\x) der = ^^"T^ — h coth w - 2 



^7r)2r $2 \^ sinh a; / sinhu; 

(43) 

To end this section, let us write down also the explicit expressions for the effective 
Lagrangians in the case of a space-time-inhomogeneous magnetic field, which is 
assumed to be small in comparison with the fermion mass scale. As follows from 
our formulae (^, (^), (|40|) and (|42|) , after performing the integration over the 
proper-time parameter, the corresponding expressions for two-derivative terms 
are: 

'''' mvrW 21 UV 7' ^ ' 

^"^^ ~ QGOttW 21 UV ^■■■)- 

These results for four dimensional QED can be compared with those in the case 
of three dimensional QED in the presence of a space-inhomogeneous magnetic 
field. For QED3 (under the same conditions) it was found 

^.er diBdiB feB\^ 
C^os _ d^m^feB^ 

These are precisely the expressions that led the authors of the paper ^ to the 
conclusion that the appearance of inhomogeneous magnetic fields is accompanied 
by a decrease of the vacuum energy. The latter would imply a vacuum instability 
which tends to generate inhomogeneous magnetic field in the vacuum state. How- 
ever, the question of vacuum stability is inaccessible in a derivative expansion, 
and in fact the infinite sum of terms with derivatives drives a system towards a 



system with uniform magnetic field p9 |. 



In the case of QED4, a space-inhomogeneous magnetic field may lead to decreas- 
ing the vacuum energy (recall that we use Minkowski metric) as well. However, 
as in the three dimensional case, it may be only an artifact of the derivative 
expansion approximation. 
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3 Conclusions 



The main result of our paper is the derivative expansion of the effective action 
for both the fermion (see formulae and (|55|)) and the scalar (see formulae (^) 
and (|4^)) electodynamics in 3 + 1 dimensions. Note that our results, in contrast 
to those obtained in the paper [Q, are derived in an explicit form for the most 
general electromagnetic field configuration. 

As for the application of the derivative expansion derived in the paper, it could be 
useful in extending the dynamical chiral symmetry breaking in QED4 P, by a 
magnetic field to the case where inhomogeneities are present^. This would allow 



contact with the GSI heavy-ion experiments |^ where strong inhomogeneous 
electromagnetic fields play an important role. 

Another application could be connected with the existence of very strong, inho- 
mogeneous fields during the electroweak phase transition in the early Universe 
1 3^ which might essentially change the character of that phase transition. 



In conclusion we note that the technique used in the present paper can be applied 
to generate the derivative expansion around covariantly constant non-abelian 
and gravitational fields. The first term in such an expansion has already been 
calculated by Savvidi ||3^ with a result generalizing Schwinger's old result [0] to 
the non-abelian gauge-field. Recently, such a calculation was extended to the 
case of nonvanishing, covariantly constant, Riemann curvature by Avramidi . 
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5 Appendix 

Here we give the functions that were used in (|35D: 
C^{a,(3) = r2tanh(ar)/7(/3r), 

C^(a,/3) = aT'H{aT)H{(3T)-— -[H{f3T) ~ H{aT)], 

— 

C|^^(a,/?,7) = — tanh(ar)tanh(/?r)i/(7r), 
8 

Cf '"(a. A ,) = ^|tanh(„.) + ta„h(^r)] ( ^("^ + W - g W _ /f^] , 

4 y a + p — 7 a + p J 

Cr^(a,/?,7) = -^tanh(ar) (^/?rg(/jr)ff(7r) - ) , 



7) = ^H{aT)HiPT)HM - ^^^^ mpr) - HM) 
^'^''^^^\Hiar)-HM)-^- ^^^^^ 



2(a + 7) ' ' ' " 2(a + 7)(a + /3) 
r iJ(/5r) r iJ(7r) 

2(a + /5)(/5-7) ^ 2(a + 7)(/3-7)' 
^ T'af3HiaT)Hif3T) r^[2{a - P) + ^ ] 



2(a-/3)(a-/3 + 7) 2(a - - /3 + 7) 

X H{-fT)[(3H{/3T)-aH{aT)] 

ar f 2(/3 + 7) 2a - /5 + 7 

+ — H[aTj 



2 ' ' V(«'-/^')(«'-7') (a-/3)2(a + 7)(a-/? + 7) 



2 ' V(a2-^2)(,32-Yi) (a_;3)2(;3_^)(a_^ + ^) 



2 'V(a^-T')(3'-7=) (a + 7)(/?-7)(a-/3 + 7), 

r 

^ 2(a-/3)(a-/5 + 7)' 



13 



-M{(3t) - — -:,H[^t) , 



+ 



(a2-^2)(y _^2)^ 
r^a fa 1 



T[//(/3r) - (7t)] / /37(/3 + 7) _ . 2\ ^(/^r 
/32-72 +2Vq; + /3 + 7 "j(Q; + /3 



)//(7t) 



/3)(a + 7) 



2(q; + /3)(q; + 7)(q; + /3 + 7) 2 1^ (a + /3) (a + 7) 



+ 



« + /5)(/5-7) (a; + 7)(/9-7) (a + 7)'(« + + 7) 
/5i7(/?r) 1 1 



(a + /3)2(a + /3 + 7) (a + + /3 + 7) (a + 7)(a + /3 + 7)7 ' 
Here we used the notation 

- xcoth(x) — 1 

= 2 ■ 
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